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We present the frequency analysis of the multicolour time series photometry of the field RRd variable BS Comae. The large 
number of data points in each of the BV(RI) C bands and the ~ 0.01 magnitude accuracy of the individual measurements 
allow us a high precision analysis of the properties of the combination frequencies due to nonlinear coupling. Through 
the combination of the frequency spectra in different colors we show that except for the components corresponding to 
the linear combinations of the two pulsation modes, there are no other components present above the millimagnitude 
amplitude level. 
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1 Introduction 

Most of the RR Lyrae variables are known to pulsate either 
in the fundamental or in the first overtone mode (RRab and 
RRc stars, respectively). There are also a few evidences for 
the existence of second mode pulsators (RRe stars) in glob- 
ular clusters (e.g. Walker & Nemec 1996) and also in the 
Galactic field (Kiss et al. 1999) but these are still considered 
to be controversial (Smith 2006). At a moderate fraction of 
RR Lyrae stars multiperiodicity can be observed. One of 
these is the phenomenon of periodic light curve modula- 
tion known as the Blazhko effect, which is known already 
for a century (see e.g. Jurcsik et al. 2005). In the case of 
the other type (denoted as RRd), simultaneous excitation of 
the fundamental and first overtone radial modes take place. 
In most cases the two types can be clearly distinguished 
observationally, since the period ratio of the latter type is 
well defined, and can vary only in a narrow interval around 
- 0.745. 

Double-mode RR Lyrae stars are very important objects 
in astrophysics for various reasons. In general, studying stel- 
lar oscillations is a good way to gain access to their phys- 
ical properties and studying RRd's is perhaps the cleanest 
field of astroseismology because unambigous mode identi- 
fication is possible only for radial pulsation. In the case of 
RRd's we have a unique opportunity to estimate absolute 
physical parameters such as the stellar mass. The schematic 
physics behind this is that radial oscillations mainly depend 
on the stellar mass and radius, and having two periods these 
parameters can be determined. Furthermore, they are also 
important test objects of nonlinear stellar pulsation models 
in the case of two radial modes (see Kollath et al. 2002). 

The first star identified as a double-mode RR Lyrae was 
AQ Leo by Jerzykiewicz & Wenzel (1977). Since then a 
large number of RRd stars were discovered mostly in Galac- 
tic globular clusters and extragalactic objects (Magellanic 
Clouds and dwarf spheroidal galaxies). Unfortunately, only 



two dozens of them are known in the Galactic field - a vast 
majority of these were found very recently by the analysis of 
the databases supplied by the automated surveys AS AS and 
NSVS (Wils 2006 and references therein). Although field 
RRd stars are of particular importance, yet most of them 
are very poorly studied. Because they are relatively bright, 
their spectroscopic follow-up study is possible with middle- 
sized telescopes. Additionally, since field RR Lyrae stars are 
distributed in a much wider metallicity interval than those 
found in clusters, they could be used to test the mass - 
metallicity distribution on a much wider metallicity range. 

2 The data 

BS Com was first considered by us as a possible target 
for our ongoing photometric survey of Blazhko stars (see 
Sodor in these proceedings) based on the investigations of 
Clementini et al. (2005). Therefore we put it on our ob- 
serving schedule. In the meantime, its double-mode nature 
was published by Wils (2006), who made rough estimates to 
the periods by the analysis of ASAS-3 data. We carried out 
four-color photometric observations with the 60 cm auto- 
matic telescope of the Konkoly Observatory, equipped with 
a Wright 750 x 1100 CCD, in the Johnson-Kron-Cousins 
BV(RI)c passbands. During the springs of 2005 and 2006 
a total number of 35 nights of observation were made. Data 
reduction and relative aperture photometry were performed 
using stanard IRAFQ packages. Our large number of data 
points (almost 1000 in V{RI) C and nearly 600 in B) with 
good phase coverage constitutes the most complete RRd 
photometric data set to date, and enable us to perform a 
Fourier-analysis of high signal-to-noise ratio. 



1 IRAF is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in As- 
tronomy, Inc., under cooperative agreement with the National Science 
Foundation. 
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Figure 1 Spectral window (upper panel) and amplitude 
spectrum (lower panel) of the V data of BS Com. The peaks 
corresponding to the radial fundamental and first overtone 
modes are marked with arrows. 

3 Frequency analysis 
3.1 The V light curve 

Since the V data have the highest accuracy among the dif- 
ferent colors, we performed a standard Fourier frequency 
analysis (ie., discrete Fourier transformation, DFT, see e.g. 
Deeming 1975) on the time series of this color, using the 
utilities of the program package MuFrAn (Kollath 1990). 
Figure Q] shows the amplitude spectrum and the spectral 
window. The fundamental and the first overtone modes and 
the ±1, 2, ... days aliases can be clearly distinguished. First 
approximation of the frequencies was obtained by a simul- 
taneous nonlinear fit of the two modes and their harmonics. 
In order to find more frequency components in the spec- 
trum, the contribution of these components were subtracted 
from the data. Folded V light curves are shown in Fig- 
ure |2] Strong variations in the residual light curves clearly 
show that the frequencies of the two pulsation modes do 
not describe the light variation within the observational ac- 
curacy. Coupling frequencies are invoked in the pulsation 
due to strong nonlinear processes in the stellar atmosphere. 
In order to identify the remaining frequency components, 
a successive prewhitening procedure was performed. The 
subsequently identified frequencies were fitted simultane- 
ously with all the previously found ones and their overall 
contribution was subtracted from the original data set. This 
iterative method increases the signal detection probability 
by eliminating the strong alias structure in the vicinity of 
a frequency peak, enabling one to identify "hidden" low- 
amplitude components suppressed by aliases and separate 
frequency peaks that are close to each other. The residual 
spectra after the subtraction of the two radial modes and 
their harmonics is shown in Figure [3^: high amplitude lin- 
ear combination components are present. The subsequent 
steps of the procedure can be followed in Figures[3]>f. As a 
result, 15 frequencies of the light variation were found and 
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Figure 2 V light curve of BS Com folded with the periods 
of the fundamental and the first overtone mode (top pan- 
els). Synthetic light curves computed from the amplitudes 
and phases of the two radial modes and their harmonics are 
shown. Middle panels show the prewhitened fundamental 
and first overtone mode light curves. Strong variations re- 
main showing the existence of additional high amplitude 
frequency components. In the bottom panels residual light 
curves after prewithening the contibution of both modes are 
shown. The presence of other strong frequency components 
is evident, which are identified as the numerous linear com- 
bination frequencies in the Fourier spectrum. Their phase 
relations are clearly denoted: the residual light variation has 
the highest amplitude when both of the modes are around 
maximum brightness. 

all but one of them were unambigously identified. 

The final set of frequencies were determined by com- 
puting numerous least squares solutions in the vicinity of 
the approximate values of /o and /i, by fixing the har- 
monics (2/o,3/o,2/i,3/i,4/i) and the 8 linear combina- 
tion frequencies shown in Figure [3] at their exact values 
and varying only the main frequencies. The solution yield- 
ing the smallest residual scatter was accepted. The full 
light curve and the corresponding Fourier fit are shown 
in Figure |U The final frequencies are /o = 2.04955 and 
/1 = 2.75432. The asymptotic error estimates of the fre- 
quencies are 0.000025d^ 1 and 0.000020a- 1 respectively. 
The 15 frequency components fit the BV(RI)c data with 
0.0119, 0.0103, 0.0099, 0.0095 magnitude rms scatter, re- 
spectively, which are near to the observational accuracy. 

3.2 Spectrum averaging 

In order to utilize the information given in the different col- 
ors, we employed the spectrum averaging method (SAM) 
as implemented by Nagy & Kovacs (2006). Considering the 
spectrum noise as random numbers, for optimizing the noise 
suppressing, the spectra were weighted inversely with their 
variances: 

* — — 2 — ; — 2 — ; — 2 — 2 — 2 a v "i — 2 — 2 a R h — 5 — T a i 
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Figure 3 Residual amplitude spectra of the V data set 
after the subsequent steps of the deconvolution procedure. 
Arrows show the identified amplitude peaks of the various 
linear combination frequency components, a: Residual af- 
ter the contribution of the modes and their harmonics were 
prewhitened from the data, b-f: The subsequent residuals 
and frequency identifications, after the removal of the pre- 
viously identified higher amplitude signals. 



where ax is the amplitude spectrum of the time series in 
the X band and ax is the standard deviation around a low- 
order polynomial fit to the red noise of the spectra. We could 
made use of SAM because we have similar numbers of data 
points in each of the colors and the photometric sampling 
was also the same that had led to very similar spectral win- 
dows. If there is any signal remaining hidden in the residual 
with not too drastically different significance in each color, 
then its detection probability increases due to the decrease 
of the noise level. Figure [5] shows the effect of averaging 
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Figure 4 V light curve of BS Com (dots) fitted with the 
Fourier solution (continous line). 



the residual spectra after all but the two smallest amplitude 
frequency components were subtracted from the time se- 
ries. The signal at 3/i — fo slightly enhanced but there is 
no effect on fo + 4/i due to its very small amplitudes in 
(RI)c- Figure|6]shows the result of the method on the resid- 
ual spectra after all the 15 frequency components were sub- 
tracted from the data. There is no additional signal over the 
millimagnitude level. The standard deviations of the SAM 
spectrum has became lower by a factor of ~ 1.5 than the 
spectra of the individual colors. This decrease is lower than 
V3, as expected from the averaging of independent random 
noise, which shows that a correlated low-frequency noise is 
present in the spectra, due to various systematics yet unfil- 
tered from the data. 



3.3 Behaviour of the Fourier-parameters 

The amplitudes of all the detected frequencies increase to- 
wards shorter wavelengths and the amplitudes of the linear 
combination frequency components are smaller than those 
of their compounding frequencies. All of the detected cou- 
pling terms are at positive linear combination frequencies of 
the modes except for f — /j and 3/i — fo- This seems to 
be in agreement with model results of Antonello & Aikawa 
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Figure 5 Residual spectra after all frequencies except 
fo + 4/i and 3/i — fo were prewhitened from the time 
series. The upper three panels show the individual spectra 
for the V(RI)c data, the bottom panel shows their average 
on a normalized scale. Whereas there is no signal increase 
°f fo + 4/i due to its very small amplitudes in (RI)c, the 
signal of 3/i — fo seems to be enhanced. 




52 



average 

jL.rifc*-ii. " j 



frequency [c/d] 

Figure 6 Residual spectra after all the 15 frequencies 
were prewhitened from the time series. The upper three pan- 
els show amplitude spectra for the V(RI)c data, the bottom 
panel shows their trend-subtracted average on the same am- 
plitude scale. The noise is suppressed by a factor of ~ 1.5 
and there is no sign of additional frequency components. A 
low-order polynomial fit to the red noise is shown in each 
panel. 



(1998). According to their nonlinear transient double-mode 
Cepheid models, frequency terms of positive linear combi- 
nations have always greater amplitudes than negative linear 
combination ones. 

Generalized amplitude ratios are defined by 

Rik = (|i| + \k\)A ik /(\i\A 10 + \k\Aoi) 

where Aik is the amplitude of the ifo + kfi component (Jur- 
csik et al. 2006). They decrease exponentially with increas- 
ing order + (Figure^), which is a similar finding as in 
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Figure 7 a: Exponential order dependence of the gener- 
alized amplitude ratios for the V data set. b: Generalized V 
phase differences plotted against the order. The relations are 
very similar for the other colors. 

the case of the triple-mode V823 Cas (Jurcsik et al. 2006). 
The generalized phase differences 



fc$oi 



where $^ is the phase of the ifo + kfi component (An- 
tonello 1994) describe phase relations independently from 
the epoch. They are quasi-coherent at a given order and they 
show linear progression with order (Figure|7j3). 

4 Summary 

We presented preliminary results of the frequency analysis 
of the double-mode RR Lyrae pulsator BS Com based on 
a substantial set of multicolor time-series photometry. In 
addition to the standard single-color Fourier analysis, we 
combined the frequency spectra in three different colors in 
order to decrease the noise level. The moderate success of 
the spectrum averaging and the presence of low-frequency 
noise indicate the importance of instrumental drift and var- 
ious systematics, yet to be filtered out from the time series. 
The full photometric time series can be fitted with 15 fre- 
quencies within the observational accuracy. Detailed values 
of the derived Fourier parameters together with the pho- 
tometric data and the estimated physical parameters from 
model results will be published elsewhere. 
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